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Abstract 

Axion-Like Particles (ALPs) are predicted by many extensions of the Standard 
Model and give rise to characteristic dimming and polarization effects in a light 
beam travelling in a magnetic field. In this Letter, we demonstrate that photon- 
ALP mixing in cosmic magnetic fields produces an observable distortion in the 
energy spectra of distant gamma-ray sources (like AGN) for ranges of the ALP 
parameters allowed by all available constraints. The resulting effect is expected to 
show up in the energy band 100 MeV — 100 GeV, and so it can be serched with the 
upcoming GLAST mission. 
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1 Introduction 



It is generally taken for granted that observations yield fair images of astro- 
nomical sources, provided sufficient care is exercised. However, environmental 
effects on the photon beam - from the source on its way to us - can mislead 
the observer, because unexpected effects can be at work. This happens e.g. 
when dust extinction and reddening become substantial, or when background 
magnetic fields affect the polarization state of radiation propagating in a cold 
plasma, thus producing a Faraday rotation. 

Remarkably enough, magnetic fields can also give rise to more subtle - and 
physically much more interesting - dimming and polarization effects in a light 
beam if photons couple to new hypothetical very light particles, to be referred 
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Fig. 1. Pair-production photon mean free path as a function of beam photon energy 
(from ref. [1]). 

to as Axion-Like Particles (ALPs). Turning the argument around, detection 
of nontrivial effects of this sort can be interpreted as observational evidence in 
favor of an ALP, thereby yielding a crucial piece of information to go beyond 
the Standard Model. 

Our aim is to show that photon- ALP mixing in cosmic magnetic fields can in- 
deed lead to the detection of ALPs in gamma-ray astronomy. More specifically, 
wc will demonstrate that photon- ALP mixing produces an observable distor- 
tion in the energy spectra of gamma-ray sources at cosmological distances - 
typically Active Galactic Nuclei (AGN) - for ranges of the ALP parameters 
which are allowed by all available constraints. The resulting effect is expected 
to show up in the energy band 100 MeV — 100 GeV, and so it can be serched 
with the upcoming GLAST mission. 

As far as the scope of this Letter is concerned, a further specification is in or- 
der. It is well known that electron-positron pair production in the scattering 
of beam photons off extragalactic background hght (EBL) becomes an impor- 
tant source of opacity whenever the corresponding photon mean free path 
is smaller than the source distance D. The energy-dependence of A-^ can be 
computed within realistic models for EBL and is reported e.g. in Fig. 1 (from 
ref. [1]). Manifestly, the resulting dimming complicates the distortion pattern 



2 



arising from photon-ALP mixing alone. In order to achieve a better under- 
standing of the latter mechanism, we find it convenient to presently discard 
EBL-induced absorption effects, deferring their analysis to a separate publi- 
cation [2]. We will therefore focus throughout on the regime in which Xj > D. 
A glance at Fig. 1 shows that this situation occurs either for beam-photon en- 
ergy E < 10^ GeV and arbitrary values of D, or else for E > 10^ GeV provided 
that the condition X-yiE) > D is explicitly enforced. 

This Letter is structured as follows. Sect. 2 offers a brief overview of the 
properties of ALPs which are of direct relevance for the subsequent discussion. 
Sect. 3 summarizes those features of cosmic magnetic fields in which observable 
photon-ALP conversion is likely to take place. Sect. 4 contains the quantitative 
estimate of the distortion of the energy spectra of AGN arising from photon- 
ALP mixing. Finally, we summarize our main conclusions in Sect. 5, where we 
also compare the proposed mechanism with similar ones recently appeared in 
the literature. 



2 Photon-ALP Mixing 



The possibility that photon mixing with a light particle alters the physical 
state of a beam was first recognized in connection with the axion, the pseudo- 
Goldstone boson associated with the Peccei-Quinn t/(l)pQ global symmetry 
invented to solve the "strong CP-problem" in a natural way [3] . In all viable 
axion models [4], the axion mass is given by m ~ 0.6 (10^GeV//a) eV, with 
fa denoting the scale at which the f/(l)pQ symmetry is spontaneously broken. 
The quark-axion Yukawa couplings induce a photon-axion interaction at one- 
loop (arising from the triangle graph with internal fermion lines), which is 
described by the effective lagrangian 

C^^-^F^'^F,^cl>^^E.Bcf>, (1) 



where (f) stands for the axion field, M~ 1.2 ■ 10^° A; (/^/lO^ GeV) GeV and 
/c ~ 1 is a parameter whose exact value depends on the specific axion model [5] 
(M is actually independent of the mass of the fermion running in the loop). 
Hence, the axion enjoys the characteristic mass-coupling relation 

^^0.7.k[l^)eY. (2) 



We stress the fact that the lagrangian C^^ naturally arises in a much broader 
class of realistic models, encompassing four-dimensional extensions of the 
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Standard Model [6], compactified Kaluza-Klein theories [7] and superstring 
theories [8]. Accordingly, C^^ is thought to describe ALPs, similar in nature 
to the axion but with m and M treated as independent parameters 0- 

A straightforward implication of C^f,^ is that the interaction eigenstates differ 
from the propagation eigenstates in the presence of a magnetic field B, thus 
generating photon-ALP interconversion; the form of C^j,^ entails that only pho- 
tons polarized in the plane containing B and the propagation direction mix 
with ALPs. As a result, a photon beam traveling in a magnetic field undergoes 
specific effects. Exchange of virtual ALPs affects the polarization state in a 
selective manner, whereas production of real ALPs - occurring for photon en- 
ergies E > m - decreases the beam intensity (besides rotating the polarization 
vector) [9]. 

Several laboratory as well as astrophysical consequences of the photon-ALP 
mixing have been addressed, in the hope to detect ALPs [10]. In particular, 
the failure to observe ALPs coming from the Sun in the CAST experiment at 
CERN has set the stark lower bound M > 1.14 ■ 10^° GeV for m < 0.02 eV [11], 
which practically coincides with the theoretical bound derived from consid- 
eration of globular cluster stars [10]. A stronger bound holds for ALPs with 
m < 10~^° eV: observations of time-lag between opposite-polarization modes 
in pulsar radio emission [12] as well as the energetics of supernova 1987a [13] 
yield M > 3 • 10^^ GeV. 

We recall that coherent photon-ALP mixing can be regarded as an oscillation 
process - much in the same way as it takes place for massive neutrinos - apart 
from the fact that an external B field is needed here, due to the spin mismatch. 

Suppose for the moment that B is homogeneous and let us denote by 
its component transverse to the propagation direction of a monochromatic 
photon beam with energy E. Then the probability that a photon will convert 
to an ALP after a distance x reads [9] 

= sin^2^ sin^ (^) , (3) 



where the photon-ALP mixing angle 6 is 
1 f \ 

= - arcsin -— — — I (4) 
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ALPs are supposed to be light enough, and for definiteness one assumes m < 1 eV. 
At variance with the axion case, here the existence of is just regarded as the 
defining feature of ALPs without bothering about its origin. 
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Fig. 2. Left panel: values of the pair (m, M) which determine the critical energy 
= I GeV, 10 GeV, 100 GeV and 1 TeV (from left to right) for a magnetic field 
strength oi B = 1 ■ 10~^ G (solid line) and B = 5 ■ 10~^ G (dotted line) and a 
plasma frequency Upi ~ 10^^'* eV. The gray region represents the values excluded 
by astrophysical arguments and by the CAST experiment. 

Right panel: same as left panel, but with B = 1-10~^ G (solid line) and B = 4-10^^ G 
(dotted line) and a plasma frequency tOpi ~ 10~^^ eV. 

and the oscillation wavenumber reads 
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so that the oscillation length is Lose = 27r/Aosc IZI- In fact, eq. (5) pertains 
to the situation in which the beam propagates in a magnetized cold plasma, 
which gives rise to an effective photon mass set by the plasma frequency 
Lj-pi = ^ATcane/fTie — 3.69 • lO^^"*^ ^ne/cm~^ eV, where Ue is the electron den- 
sity {rrie denotes the electron mass). 

A deeper insight into the physics of photon- ALP oscillations can be gained by 
introducing the critical energy 



\m 



2Bi 



0.26-10 



10 



uo. 



;i0-i0eV)2 \ Bt ) VlOiOGeV 



M 



eV (6) 



(see Fig. 2). For further needs, we note that eqs. (3) and (5) can be rewritten 



2 Since we are dealing with weak magnetic fields, their contribution to the vacuum 
refractive index is negligible [14]. 
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and 
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respectively. Accordingly, the situation can be schematized as follows. 

Strong-mixing regime ~ In the high-energy limit E ^ E^, we immediately have 
Aosc — Bt/M, the photon- ALP mixing is maximal {9 ~ 7r/4) and the conver- 
sion probability becomes energy-independent. 

Weak-mixing regime - In the opposite low-energy limit E ^ E^, we get Aqsc — 
~ — LJpi\/2E. The mixing is small {6 <^ 1), photon-ALP oscillations be- 
come dispersive - since now both the mixing angle and the oscillation length 
are energy-dependent - and their amphtude gets reduced by the factor [E/E^y. 

In cither case, the simpler behavior P^\^[x) ~ {Bt x/2My emerges for an 
oscillation length Lqsc ^ x. 



3 Cosmic magnetic fields 



As already stressed, photon-ALP conversion requires the presence of a mag- 
netic field playing the role of a catalyst. Below, we consider those cosmic 
magnetic fields which are likely to affect in a substantial manner the physical 
state of a photon beam from a distant AGN (within the ALP scenario outlined 
in Sect. 2). 

Generally speaking, the origin and structure of magnetic fields in the Uni- 
verse is still unknown. A possibility is that very small magnetic fields present 
in the early Universe were subsequently amplified by the process of structure 
formation [15]. An alternative option is that magnetic fields have been gener- 
ated in the low-redshift Universe by energetic quasar outflows [16]. Finally, it 
has been suggested that the seeds of extragalactic magnetic fields originated 
from the so-called Biermann battery effect [17], namely from electric currents 
driven by merger shocks during structure formation processes. Presumably, 
all these effects can take place, even if it is presently impossible to establish 
their relative importance. 
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Observations show that cosmic magnetic fields have a complicated morphol- 
ogy, which evidently reflects both the pattern of baryonic structure formation 
and its subsequent evolutionary history [18]. In spite of the fact that they 
come in a wide variety of configurations and strengths, many of them cannot 
be considered as uniform over the typical distances traveled by photons from 
cosmological sources. As a consequence, one cannot evaluate the photon- ALP 
transition probability by blind application of eq. (3). 

It turns out that - at least to a first approximation - one can assume that 
nonuniform cosmic magnetic fields B have a cellular structure (more about 
this, later). That is, B is supposed to be constant over a domain of size Ldom 
equal to its coherence length, with B randomly changing its direction from 
one domain to another but keeping approximately the same strength. Over 
distances D ^ I/dom, the actual conversion probability P^^^{D) arises as the 

incoherent average of ^^^^(-Z^dom) over the N ~ {D/L^om) domains crossed by 
the beam. One finds [19] 



(9) 



which can be approximated as P^^^{D) ~ 0.5 ■ A^P^j|.0(Ldom) for N and Ldom 
such that N P^^_l^{Ldom) 1- Alternatively, P^^^{D) saturates in the limit 
N Pl'^_ls{L^om} ^ 1, SO that on average one-third of the photons become ALPs. 



3. 1 Large-scale magnetic fields 



So far, observations have failed to detect the existence of magnetic fields over 
cosmological scales and only upper limits are available on their strength and 
coherence length. Typically, one gets B < 10~^ — 10"^ G over Megaparsec 
scales [18,20]. 

In the lack of any reliable information, we will carry out our analysis for 
B — 1 ■ 10~^ G and B — 5 ■ 10~^ G, and -Ldom — 1 Mpc, which are close to 
existing upper limits but consistent with them. 

Still, it is interesting to notice that our preferred values are suggested by 
a simple heuristic argument [21]. Observations yield B > 10~^G in col- 
lapsed baryonic structures with overdensity 5 ~ 10^. Flux conservation during 
gravitational collapse (adiabatic compression) entails B ~ 5^/^. So, we get 
B > 10~^ G in the intergalactic medium. Moreover, in the quasar outflow 
model the cellular structure of the magnetic fields emerges naturally, with a 
coherence length of a Megaparsec scale. 
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Besides from magnetic fields, the physical state of the beam propagating over 
cosmological distances is also affected by the presence of a cold plasma in 
intergalactic space. The absence of the Gunn-Peterson effect is usually taken as 
an evidence that the intergalactic medium is ionized with rig ^ 10~^cm~^ [22], 
resulting in the plasma frequency cupi ~ 1.17 • 10~^^ eV, in agreement with the 
WMAP upper bound Ue < 2.7 • 10~^cm~^ on the baryon density [23]. 

3.2 Intracluster magnetic fields 

A better situation concerns clusters of galaxies. Indeed, observations have 
shown that the presence of magnetic fields with average strength B ~ 10^^ G 
is a typical feature of the intracluster region. Somewhat stronger values are 
detected in the cores of regular clusters. Even more remarkable is the fact 
that observations are able to yield information about the associated coher- 
ence length, which turns out to be of the order of 10 kpc [24]. A cellular 
structure for the intracluster magnetic field is usually assumed, with domain 
size I/dom — 10 kpc. 

Just as in the previous case, plasma effects are expected to show up when 
the beam crosses a cluster. Specifically, the electron density of the intra- 
cluster medium is — 1.0 • 10~^cm~^ [25], which yields a plasma frequency 
a;pi~ 1.2- 10-12 eV. 

3.3 Galactic magnetic fields 

Observations over the last three decades have led to a rather detailed picture 
of the magnetic field in the Milky Way. Perhaps, the most important feature 
of the Galactic magnetic field is that it consists of two components. 

Regular component - Measurements of Faraday rotation based on pulsar ob- 
servations have shown that this component is parallel to the Galactic plane. 
Its strength varies between i? ~ 2 ■ 10"^ G in the Solar neighbourhood and 
B ~ 4-10-^ G at 3 kpc from the centre [26]. Moreover, the associated coherence 
length is of the order of 10 kpc. 

Turbulent component - Over much smaller scales, the dominant Galactic mag- 
netic field appears to be stochastic, with a Kolmogorov spectrum a — 5/3 [27]. 
In practice, this component can be described by a cellular structure, with 
strength S ~ 1 • 10"^ G and domain size -Ldom — lO^^ pc. 

Inside the Milky Way disk the electron density is rig — 1.1 • 10-^ cm"^ [28], 
which gives a plasma frequency cUpi ~ 4.1 • 10"^^ eV. 
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4 Spectral distortion 



Our proposal concerns the distortion of the energy spectra of extragalactic 
sources like AGN as induced by photon- ALP conversion in intervening mag- 
netic fields (of the kind discussed in Sect. 3). Presently, the source emission 
spectrum dN/dE gets modified along the line-of-sight in such a way that 
at the observer position it becomes dN/dE times the total photon survival 
probability P^^^{D). Because P^^^{D) = 1 — P^^^{D), we see that the ob- 
served spectraA distortion is just dN / dE ■ P^^^(D) . That is to say, the emission 
spectrum merely gets distorted in proportion to the photon-ALP conversion 
probability regardless of the actual spectral shape. We stress that this circum- 
stance greatly simplifies our analysis, since it dispenses us from committing 
ourself with a specific source spectrum. 

Owing to eq. (9), the size of the observed spectral distortion increases with 
D. Yet, a larger D both makes the source fainter and enhances the EBL- 
induced absorption at high energies. In the analysis to follow, we will adopt 
for definiteness the realistic values D = 200 Mpc and D = 500 Mpc whenever 
necessary. 

It goes without saying that the spectral energy distortion has to be mea- 
sured well enough in order to disentangle the effect in question from other 
uncertainties. It should also be kept in mind that a large part of the error in 
very-high-energy gamma-ray detectors is correlated [29], so that the ratio be- 
tween the yields in two different energy points can be measured with a relative 
uncertainty of order 10% in a Imaging Atmosphreric Cherenkov Telescope and 
well below 10% in GLAST. Thus, a spectral distortion larger than 10% will 
be regarded as observable throughout the subsequent discussion. 

We proceed to investigate the behavior of the photon-ALP conversion prob- 
ability as a function of Over a single magnetic domain, P^^^(Ldom) is 
computed from eq. (7) and is shown in the left panels of Figs. 4 and 6 for 
suitable values of the magnetic field. Over the whole distance D = N Ldom, 
P^^^{D) is given by eq. (9) and is similarly illustrated in Fig. 5 and in the right 
panels of Figs. 4 and 6. We see that the conversion efficiency increases with 
energy as long as E < E^, while becomes maximal for E > E^^. A character- 
istic feature shows up due to the drastic change in the yield when comparing 
energies above and below E^ by two-three orders of magnitude - this is indeed 
the signature of the effect we are looking for 0. 

We identify the energy band between 100 MeV and 100 GeV as the best 
compromise between the detector sensitivity and the lack of EBL-induced 

^ In agreement with the discussion in Sect. 2, an oscillatory pattern is present for 
E < E^, until it becomes unobservable at sufficiently low energy. 
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Fig. 3. Left panels: region in the (m, M) space which determine the critical energy 
between 100 GeV and 1 TeV for a magnetic field strength of B = 1 ■ 10~^ G (upper 
plot) and B = 5 ■ 10~^ G (lower plot) and a plasma frequency Wpi ~ 10~^^ cV. The 
dark gray region represents the values excluded by astrophysical arguments and by 
the CAST experiment. 

Right panels: same as left panels with however B = 1 ■ 10~^ G (upper plot) and 
B = 4 - 10~^ G (lower plot) and a plasma frequency LOp\ ~ 10~^^ eV. 

absorption. It is easy to check that condition \{E) > D is presently met. 
The energy band in question is almost completely unexplored at present, but 
it will soon become accessible with the GLAST sateUite. 

According to the foregoing discussion, observable effects can be detected in 
such an energy band provided the critical energy i?* lies just above its upper 
edge, namely for ~ 10^ GeV — 1 TeV. The constraints implied by the latter 
condition on the parameters m and M are reported in Fig. 3 for suitable values 
of the magnetic field. Correspondingly, we find that in the allowed region of the 
(m, M) space the less stringent CAST bound on M applies and that plasma 
effects are unimportant. 

Below, we evaluate the spectral energy distortion for a distant AGN as pro- 
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Fig. 4. Conversion probability versus photon energy in units of E^, in the large-scale 
magnetic field. The plots in the left panels show the conversion probability over 
a single magnetic domain p!^\^, whereas those in the right panels represent the 
total conversion probability P-y^ff, over N = 200 magnetic domains. The plots in 
the upper panels are obtained for M = 3 • lO^'^ GeV, while those in the lower panels 
arise for M = 1 • 10^^ GeV. Dotted and solid lines correspond to -B = 1 ■ 10~^ G 
and B = 5 ■ 10~^ G, respectively. 

duced by photon- ALP conversion occurring in the three magnetic environ- 
ments considered in Sect. 3. Clearly, all we have to do is to evaluate P^^^{D). 
Observable effects get singled out by the requirement P^^^{D) > 0.1. 



Large-scale contribution 



The conversion probabihty in the large-scale magnetic field is shown in Fig. 4 
and Fig. 5. We have computed the effect for source distances D — 200 Mpc 
(Fig. 4) and D = 500 Mpc (Fig. 5), which correspond to = 200 and N = 500 
magnetic domains crossed by the beam, respectively. The left panels of Fig. 4 
show the effect over a single domain of the magnetic field (this is obviously 
independent of distance). We have considered two different values of the in- 
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Fig. 5. Conversion probability P-y^c/, versus photon energy in units of in the 
large-scale magnetic field: same as in the right panels of Fig. 4 but over N = 500 
magnetic domains. The left plot is obtained for M = 3 • 10^'' GeV, while the right 
one arises for M = 1 • 10^^ GeV. Dotted and solid lines correspond to S = 1 • 10~^ G 
and B = 5 ■ 10~^ G, respectively. 



verse photon coupling constant M (M = 3 • 10^° GeV and M = 1 ■ 10^^ GeV), 
and two different values of the magnetic field strength B {B = 1 ■ 10~^ G and 
5 = 5-10-'' G). 

We see that for M sufficiently close to the CAST lower bound (upper panels 
in Fig. 4 and left panel in Fig. 5) the spectral distortion is observable provided 
B is roughly within one order of magnitude from the upper limit. A similar 
result shows up in the opposite situation, namely for B sufficiently close to 
the upper bound (solid lines) and M roughly within one order of magnitude 
from the CAST lower limit (lower panels in Fig. 4 and right panel in Fig. 5). 
This conclusion is practically unaffected by the source distance. 



4-2 Intraduster contribution 



The conversion probability in the magnetic field of a galaxy cluster is shown 
in Fig. 6. We assume that the cluster has a typical size of 1 Mpc, so that 
= 100 magnetic domains are crossed by the beam. These plots correspond 
to M = 3 ■ 10^° GeV and M = 1 • 10" GeV, and to £ = 1 ■ lO'^^ G and 
S = 4 • 10-6 G. 



Clearly, if the beam goes through a cluster of galaxies the spectral energy 
distortion turns out to be observable for all our preferred values of M and 
B, that is to say provided these parameters lie roughly within one order of 
magnitude from their bounds. 



12 



1 1 I 0.35 




lO"* 10"^ 10"' 10° 10' 10^ 10'° 10"^ 10'' 10° 10' 10^ 

E/E. E/E, 



Fig. 6. Conversion probability versus photon energy in units of £'* in intracluster 
magnetic fields as well as in the regular Galactic field. The plots in the left panels 
show the conversion probability over a single magnetic domain P^^^^, whereas those 
in the right panels represent the total conversion probability P-^,^^ over N = 100 
magnetic domains. The plots in the upper panels are obtained for M = 3- 10^° GeV, 
while those in the lower panels arise for M = 1 • 10^^ GeV. Dotted and solid lines 
correspond to 5 = 1 ■ 10~^ G and B = 4 - 10~^ G, respectively. 

4-3 Galactic contribution 

The situation concerning photon- ALP conversion in the Milky Way can be 
summarized as follows. 

Regular component - The conversion probability in this magnetic field is com- 
puted directly from eq. (7) and it is shown in the left panels of Fig. 6. These 
plots again correspond to M = 3 • 10^° GeV and M = 1 • 10^^ GeV, and to 
B = 1-10~^ G and B = 4-10^^ G. The situation is analogous to what we found 
in the case of the large-scale magnetic field. That is, for M close enough to the 
CAST lower bound the spectral distortion is observable provided B is roughly 
within one order of magnitude from the upper limit. Similarly, observability 
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is ensured for B sufficiently close to the upper bound and M roughly within 
one order of magnitude from the CAST lower limit. 

Turbulent component - It is straightforward to realize that we presently have 
Lose 3> Ldom for all experimentally allowed values of m and M. Therefore, 

P^°I^(^dom) ^ (5TLdom/(2M))2 and we get P^^^ < 2.3 ■ IQ-^^ by enforcing 
M > 10^° GeV. Because here the number of magnetic domains is ~ 10®, we 
end up with P-y^^ < 10^^. Thus, we see that in this case no spectral energy 
distortion is observable. 



5 Discussion and Conclusions 

We have presented a mechanism whereby the energy spectra of gamma-ray 
sources at cosmological distances get distorted due to photon- ALP mixing 
taking place in the cosmic magnetic fields crossed by the beam on its way 
to us. We have attempted to identify the ranges of the ALP parameters for 
which the effect in question can be observed with GLAST. Unfortunately, 
the uncertainties in the properties of cosmic magnetic fields prevent us from 
making sharp statements about an exclusion plot in the ALP parameter space. 
Nevertheless, we have succeeded in showing that observability is achieved for 
ranges of the ALP parameters which are allowed by all available constraints. 

Large-scale magnetic fields as well as the regular Galactic component turn out 
to be nearly equally efficient at producing an observable distortion. Whenever 
the beam crosses a cluster of galaxy, its intracluster magnetic field is even more 
efficient in that respect. The latter circumstance suggests to look at similar 
sources however in different directions, so that cluster crossing occurs only 
for one line of sight. Directionality can also be instrumental in detecting the 
spectral enery distortion due to the Galactic regular magnetic field, since its 
morphology is presently fairly well known. 

Our proposal shares some similarities with the one advanced a few years ago 
by Gsaki, Kaloper and Terning (GKT) [30] as an explanation for the observed 
dimming of distant type la supernovae [31]. Currently, such a dimming is inter- 
preted as evidence for an accelerated cosmic expansion, presumably triggered 
by a mysterious dark energy [32]. Instead, GKT suggested that the supernovae 
under consideration look fainter than expected simply because some photons 
en route to us become ALPs in extragalactic magnetic fields, thereby escaping 
detection. Unfortunately, subsequent studies have shown that this proposal 
gets ruled out for almost all values of the parameter space [33]. In particular, 
plasma effects make the dimming of type la supernovae excessively chromatic. 
We stress that these problems are automatically avoided in our case simply 
because any dimming effect disappears at energy £■ ^ ~ 10^ GeV — 1 TeV. 
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A somewhat different idea has recently been put forward by Hooper and Ser- 
pico [34], and by Hochnmth and Sigl [35]. Here, it is proposed that photon- 
ALP conversion can take place inside gamma-ray sources, thanks to their 
strong magnetic fields. We have seen in Sect. 2 that efficient photon-ALP 
conversion requires the strong-mixing regime to be realized. Moreover, it fol- 
lows from eq. (3) that P^^^{x) becomes maximal for Aqsc^:^ ~ 1, namely for 
Bt/M X ~ 1 (this just follows from eq. (5) in the strong- mixing case). The 
latter condition is similar the Hillas criterion [36] concerning the acceleration 
of cosmic-ray particles, and this circumstance is used to argue that such a 
regime should take place in some astrophysical sources. Manifestly, both this 
mechanism and the one proposed in this Letter can be operative at the same 
time. 

We thank Sasha Dolgov and Dario Grasso for discussions. 
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